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Abstract 

CFD  examples  at  ADD  are  introduced  to  show  their  variety  at  its  application  in  the  course  of  a  missile  design.  Four 
examples  are  an  ogive-cylinder  and  boat  tail,  nose  spike,  vertical  launcher  internal  and  side  jet  interaction  flows  at 
supersonic  flow  region.  Various  means  of  validation  for  those  complex  flows  are  also  described.  This  paper  is  thus 
intended  to  show  how  CFD  and  its  validation  share  their  role  at  the  ADD  aerodynamic  research  laboratory. 
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1.  Introduction 

In  the  course  of  a  missile  system  design,  generation  of  aerodynamic  performance  data  comes  first  to  scale  the 
configuration.  For  this  activity  a  standard  empirical  or  semi-empirical  code  is  utilized  to  generate  a  quick 
performance  data.  Due  to  the  disparity  between  the  input  geometry  for  these  empirical  codes  and  the  real  flow 
conditions  or  the  geometry,  wind  tunnel  testing  is  done  to  supplement  the  aerodynamic  data  for  evaluation  of  flight 
simulation.  However,  there  are  other  situations  which  demand  a  detailed  flow  simulation  for  phenomenological 
understanding.  Cases  in  point  are  supersonic  missile  body  flow  at  high  angle  of  attack,  boat  tail  and  spike  design,  jet 
impingement  from  rocket  plume  on  an  opposing  wall,  and  side  jet  interactions  over  a  supersonic  missile.  In  this 
introductory  paper,  a  few  examples  for  CFD  applications  that  are  encountered  during  the  missile  design  are  given 
along  with  validation  methods  for  the  computed  results. 

Numerical  study  of  the  jet  flow  is  carried  out  using  a  three-dimensional  Navier-Stokes  solver  incorporating  the 
Characteristic  Flux  Difference  Splitting  (CFDS)  method  [1].  CFDS  method,  a  variant  of  Roe's  flux-difference 
splitting  formulation,  has  been  developed  formerly  as  a  viable  engineering  prediction  tool  for  aerodynamic  design 
and  had  indeed  shown  its  versatility  in  computing  complex  flows  [2].  The  CFDS  method,  however,  being  rooted 
from  Roe’s  flux-difference  is  also  susceptible  to  the  carbuncle  problem.  This  shock  instability  problem  has  been 
cured  and  the  code  is  then  applied  to  many  shock/boundary  interaction  situations  cited  herein. 

2.  Numerical  Method 

The  Characteristic  Flux  Difference  Splitting  (CFDS),  numerical  method  for  the  three-dimensional  Navier-Stokes  has 
been  applied  to  various  complex  flows  and  validated  over  the  past  few  years  [1].  The  governing  Navier-Stokes 
equations  employed  in  the  generalized  coordinate  system  ( 77,  ^ )  are  expressed  for  the  conservative  variable 
vector  as 

1  dO  d  ^  ^  d  d  ^  ^ 

at  og  or]  0(p 
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where  J  ^  is  the  Jacobian  of  the  transformation,  Q  is  the  conservative  variable  vector,  F ,  G  and  H  are  inviscid 

flux  vectors,  and,  Py’Gv  and  are  viscous  flux  vectors.  The  inviscid  fluxes  are  linearized  and  split  for  upwind 
discretizations  by 

A^F  =  AAQ  =  (A^  +  A  )AQ  and  A"  =  MfAW'T-*  (2) 


yielding 

J  '  +  A  Ajg  +  +  B  A,Q  +  C*V,2  +  C  =  0 


(3) 


where  SQ  =  and  the  overbar  means  the  associated  variable  is  space-averaged  over  the  interval,  [j, 

j+l],  M  or  M"‘  is  a  transformation  matrix  between  the  conservative  variable  vector  Q  and  the  primitive  variable 
vector,  Q  .  T  or  its  inverse  T  ^  is  defined  to  be  a  transformation  matrix  between  the  primitive  variable  vector 

Q  and  the  characteristic  variable  vector,  Q  . 

In  the  case  of  supersonic  jet  impingement  calculations,  the  entropy  fixing  formula  employed  are 

I  A  |=  )  if  \A\<£  (4) 

2^ 

with  S  =constant  to  prevent  shock  instability  problem. 


3.  CFD  Applications  in  Aerodynamic  Design  and  Analysis 
3.1  Supersonic  Ogive-cylinder 

Supersonic  flow  past  an  ogive-cylinder  with  the  secant  ogive  of  3-caliber  length  is  run  for  Mach  3  with  a  =10.0° 

and  Re=6xl0^[3].  Degani-Schiff  [4]  turbulent  model,  modified  from  Baldwin-Lomax  model,  was  used  since  the 
flow  involves  cross  flow  separation  at  high  angle  of  attack.  Pressure  distributions  on  the  upper  and  lower  surface 
along  the  longitudinal  body  are  compared  with  experimental  data  [3]  in  Fig.  1  with  a  good  match  even  with  1st- 
order  differencing.  Although  Ist-order  solution  follows  the  experimental  trend  pretty  well,  the  second-order  spatial 
accuracy  improves  the  quality  of  the  solution  as  is  apparent  from  Fig.  2  which  compares  circumferential  pressure 
distributions  at  x=5.77  caliber  location. 
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Figure  1.  Jet  Surface  pressure  distribution 

o 

over  an  ogive  at  M=3.0,  a  =10.0  . 


Figure  2.  Circumferential  pressure  distribution 
at  x=5.77D  for  M=3.0,  a  =10.0  . 


3.2  Boat  Tail  Configuration  Design  and  Nose  Spike  Design 

Boat  tail  design  of  a  missile  is  a  big  issue  concerning  the  reduction  of  drag  of  the  missile.  Drag  changes  are 
estimated  as  a  function  of  Sr,  Sr=Sbase/Snozzle.  Figure  3  shows  Mach  contours  around  the  missile  base  with  boat 
tail  angle  16  degree  with  Sr=1.45.  The  jet-off  simulation  result  is  compared  on  the  upper  half  and  the  jet-on 
simulation  result  is  plotted  on  the  lower  half  as  shown  in  the  Fig.  3.  In  case  of  jet-off,  recirculation  zone  is  appeared 
near  the  base  accompanied  by  the  shear  layer  adjacent  to  this  recirculation  zone.  A  Mach  disk  at  about  0.3D  distance 
from  the  base  can  be  seen  in  the  jet-on  case  because  this  flow  is  under  expanded. 

The  jet-off  base  drag  (blue  dashed  line  with  open  rectangle  symbol)  goes  up  and  the  jet-off  boat  tail  drag  (blue 
dashed  dot  line  with  open  rectangle  symbol)  goes  down  with  the  increase  of  Sr  as  shown  in  Fig. 4.  The  jet-on  drag 
thus  steadily  increases  when  the  area  ratio  is  increased.  Both  the  base  drag  and  boat  tail  drag  when  jet  is  on  increase 
until  the  area  ratio  reaches  1.0. 

Unsteady  turbulent  numerical  simulation  is  done  using  some  spike  configurations  in  front  of  missile  nose  in  an 
effort  to  reduce  the  drag.  Pressure  distributions  and  stream  lines  are  shown  in  Fig.  5.  Total  drag  reduction  of  30 
percent  at  low  angle  of  attack  is  gained  due  to  the  break-up  of  bow  shock  from  the  spike.  The  calculated  pressure 
oscillation  is  estimated  at  the  frequency  range  of  800-1000Hz.  The  pressure  oscillation  frequency  is  measured  from 
the  flight  test  with  the  range  of  900-1000Hz. 
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<  Mach  Number  Contours  > 


Jet  OFF 
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Figure  3.  Mach  contours  with  jet  on/off, 
Sr=1.45,  M=0.85,  boat  tail  angle=16  deg. 


Figure  4.  Drag  curves  for  varied  Sr,  M=0.85, 
hoat  tail  angle=16  deg. 


Figure  5.  Unsteady  turbulent  flow  analysis  at  M=1.2. 


3.3  VLS  Internal  Flow  Simulation  with  supersonic  jet  impingement  onto  wall 

The  VLS  design  challenge  is  to  contain  the  initial  impact  of  the  jet  plume  and  safely  discharge  the  rocket  exhaust 
gas  during  launch  of  the  missiles.  The  challenge  needs  innovative  mechanical  design  and  extensive  use  of  thermal 
protection  materials  that  require  understanding  of  the  flow  structures  of  the  exhaust  plumes.  When  the  jet  plume 
exhausts  from  rocket  motor,  the  jet  flow  impinges  on  the  bottom  wall  and  then  exits  through  the  uptake  after 
circulating  in  the  plenum  chamber. 

To  study  the  supersonic  jet  impingement  inside  the  VLS  plenum,  flat  plate  jet  impingement  test  case  is  simulated 
and  validated  experimentally  at  first.  The  computational  geometry  used  for  jet  impingement  onto  flat  plate  is  shown 
in  Fig.  6.  The  nozzle  is  mounted  perpendicular  to  the  flat  plate.  The  computational  domain  starts  from  the  nozzle 
throat  with  Mach  1.0  condition.  The  boundary  conditions  of  this  nozzle  throat  are  calculated  from  isentropic 
relations  and  perfect  gas  law.  Strong  normal  shock  is  formed  over  the  plate  when  a  supersonic  jet  plume  exhausts 
against  the  plate.  If  the  grid  system  used  in  numerical  computation  is  aligned  with  this  normal  shock,  the  shock 
instability  occurs.  Figure  7(a)  shows  Mach  number  contours  in  symmetric  plane  contaminated  with  shock  instability, 
so-called  “carbuncle  problem”.  This  shock  instability  is  cured  by  fixing  near  zero  eigenvalue  in  the  numerical 
dissipation  term,  yielding  the  stable  results  as  shown  in  Fig.  7(b). 
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(a)  without  eigenvalue  fixing  (b)  with  modified  eigenvalue 
Figure  7.  Mach  contours  in  svmmetric  plane  for  H=4D. 


Frequency(Hz) 


Figure  8.  Spectra  for  H=5.4D. 


Figure  9.  Sensor  locations  on 
the  flat  plate. 


Figure  10.  Ablation  of  flat 
plate  after  motor  test. 


Ablative  Material 
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Figure  11.  Sketch  of  motor  test  facilities  (left)  and  actual  test  scene  (right). 


It  is  observed  that  the  pressure  oscillates  with  certain  frequency.  The  pressure  fluctuation  at  the  center  position 
reveals  dominant  frequency  of  2.8kHz  through  the  Fourier  transform  procedure  as  shown  in  Fig.  8.  Pressure  and 
temperature  sensors  are  located  on  the  flat  plate  with  15  mm  space  interval  as  shown  Fig.  9  with  seven  pressure 
holes  and  four  temperature  gages.  Figure  10  reveals  the  magnitude  of  ablation  on  the  flat  plate  after  the  motor  test. 
Figure  1 1  illustrates  the  test  motor  and  the  test  setup. 

The  proven  methodology  on  the  supersonic  jet  impingement  flow  is  then  applied  to  VLS  launcher.  The 
computational  grid  system  for  flow  simulation  of  VLS  with  1 8  blocks  are  overlapped  with  neighboring  blocks  as 
shown  Fig.  12(a).  Flow  simulations  of  only  two  motor  cases  are  sufficient  because  of  geometric  similarity  even 
though  total  of  eight  missiles  are  housed  in  VLS.  Figures  12. (b)  and  12. (c)  show  Mach  contours  and  pressure 
contours  of  yz-plane  view  and  xy-plane  view,  respectively,  in  the  plane  cutting  through  the  motor. 


(a)  Overall  grid  (b)  Mach  contours  in  yz-plane  (c)  Pressure  contours  in  yz-plane 

Figure  12.  Computational  grid  topology  and  Mach  contours  for  VLS  analysis. 


Figure  13  shows  the  pressure  distributions  on  the  plenum  bottom  with  0.5  atm.  interval.  The  maximum  pressure  is 
about  20  atm.  as  illustrated  in  Fig.  13.  Since  the  primary  interest  is  whether  the  jet  plume  is  exhausted  without 
causing  too  much  erosion  on  the  bottom  as  well  as  on  the  side  wall,  actual  test  firing  was  performed  to  partially 
verify  the  flow  pattern  as  well  as  the  ablation  thickness.  Pressure  measure  points  at  the  plenum  top  plate  are  denoted 
in  Fig.  14.  The  measured  data  and  simulated  results  are  compared  in  Fig.  15  and  Fig.  16.  The  measured  peak 
pressure  at  PI  which  is  at  the  nearest  position  to  the  motor  in  Fig.  15  is  greater  than  the  simulated  pressure.  The 
pressure  at  P5  shows  good  match  between  experimental  and  CFD  data. 
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Figure  13.  Pressure  contours  of  plenum  bottom 
wall  with  AP  =0.5  (atm.). 
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Figure  15.  Pressure  PI  comparisons. 


Figure  14.  Sensor  locations  on  the  plenum  top 
plate. 


Figure  16.  Pressure  P5  comparisons. 


3.4  Aerodynamic  Modeling  of  Side  Jet  Influence  on  Body-Tail  Missile 
Configuration 

Side  Jets  are  located  on  the  forebody  section  to  cause  sudden  changes  in  flight  path.  The  location  on  the  missile 
body  is  shown  and  a  typical  jet  interaction  flow  pattern  is  summarized  in  Fig.  17.  Figure  17  contains  Mach  contours 
in  the  symmetry  plane  when  the  side  jet  is  turned  on,  capturing  a  strong  bow  shock  and  barrel  shock,  among  other 
features.  Both  the  pressure  contours  on  the  body  surface  and  the  wall  pressure  distribution  near  the  jet  exit  plane 
illustrate  severe  changes  in  the  pressure  field.  The  exit  flow  condition  of  the  side  jet  thruster  is  obtained  from  a 
separate  simulation  and  reveals  a  very  severe  pressure  condition.  The  jet  exit  pressure  is  about  500psi  and  exit  Mach 
number,  2.4.  The  pressure  ratio  of  the  exit  pressure  to  upstream  pressure  is,  for  example,  137  at  10km  altitude.  This 
ratio  will  increase  as  the  flight  altitude  goes  up  and  eventually  pose  a  difficulty  in  conducting  wind  tunnel  test. 


A  key  question  prior  to  wind  tunnel  test  was  to  uncover  the  similarity  law  for  the  side  jet  conditions  between  the 
flight  and  wind  tunnel  models.  The  similarity  parameters  are  i)  static  pressure  ratio,  ii)  momentum  flux  ratio 
between  free  stream  and  jet  exit  plane,  and  iii)  temperature  ratio.  However,  the  single  most  important  similarity 
parameter  may  be  the  ratio  between  the  upstream  pressure  and  pressure  at  the  jet  exit  plane.  Primary  aim  of  current 
side  jet  analysis  is  to  provide  6DOF  aerodynamic  module  for  guidance  and  control  simulation  with  and  without  side 
jet.  However,  the  influence  is  likely  a  function  of  flight  Mach  number,  altitude,  angle  of  attack,  orientation  of  jet,  tail 
deflection  and  the  number  of  side  jets  operated  at  one  time. 


After  further  CFD  simulations  and  comparisons  of  data  with  data  obtained  from  in-house  wind  tunnel  tests, 
empirical  formula  which  accounts  for  the  effect  of  side  jet  force  are  defined  as: 


c  =c 


N 


+  K  N 


Without  Jet 


'Ciu  ^  Jet 


T 

Single  Jet 
Qref^ref 


c  =c 


M 


T  T 

+  k  N  ^^^SleJet  J^jet 


Without  Jet  '^M 


Cm  Jet 


QrefSyef  Cef 


(5) 

(6) 
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where, 

K,  =  H,  Nj^, ,  )  (7) 

=  8  {a,  Mach,  H ,  N  )  (8) 

where  the  empirical  factor  K  was  introduced  for  force  and  moment  coefficients.  It  is  called  hereafter  the  jet 
effectiveness  factor,  K.  This  study  thus  explores  the  effect  of  side  jet  on  pressure  changes  on  the  body-tail 
configuration  or  aerodynamic  coefficients  through  the  changes  of  this  value  from  1.0.  The  aerodynamic  changes 
influenced  by  side  jet  activation  can  thus  be  quantified  through  the  K  factor  for  all  components  of  forces  and 
moments.  If  there  is  no  side  jet  influence  on  the  wall  pressure  changes  on  the  body  and  tail.  If  it  is  greater  or  less 
than  1.0,  that  difference  represents  the  amount  of  influence  by  the  side  jet. 


MACH  CONTOUR  (H=1  OKm,  Mach  2.6,  Aoa  0^ 


SURFACE  PRESSURE  DISTRIBUTION 

(Poe/Poo=137) 


PRESSURE  DISTRIBUTION  ALONG  THE  JET-CENTER  LINE 


Single  Sidejet  Thruster 

Divergence  Angle  ;  8  deg. 


Figure  17.  Jet  interaction  flowfield  near  and  on  the  body  surface. 

Having  studied  the  behavior  of  the  jet  effectiveness  factor  K  computationally,  wind  tunnel  test  is  still  performed  to 
map  out  the  overall  trend  of  the  K  factor  for  all  flow  and  geometry  parameters  since  the  CFD  alone  can  not  cover  the 
whole  flight  envelop.  After  much  planning  and  preparation,  two  sets  of  wind  tunnel  tests  have  been  conducted  at  the 
S2  and  S3  wind  tunnel  facilities  of  ONERA  at  Modane  in  France  March  2003  and  April  2007.  A  test  matrix  was 
made  to  cover  Mach,  altitude,  angle  of  attack,  jet  orientation  or  bank  angle,  and  number  of  jets.  Figure  18  shows  a 
schlieren  picture  of  flow  field  at  Mach  2.3,  pressure  ratio  of  150  with  the  jet  on  the  leeward  side.  It  shows  nose 
shock,  strong  bow  shock,  barrel  shock,  expansion  wave  at  the  nose  junction  and  a  rather  large  wake  region. 

Results  obtained  from  S2  tunnel  are  reprocessed  to  extract  the  jet  effectiveness  factor  K.  The  K  curves  are  then 
grouped  and  plotted  for  various  flow  parameters.  Each  curve  is  plotted  as  a  function  of  angle  of  attack. 

In  general,  strong  interaction  between  the  jet  and  the  body  is  observed  when  the  jet  is  activated  on  the  windward  side 
(Bank=180).  The  wall  pressure  on  the  body  downstream  of  the  jet  goes  up  since  the  body  shock  generated  from  the 
nose  confines  the  flow  within  the  shock  layer  and  thereby  keeping  the  flow  attached.  All  the  K  curves  with  the 
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windward  jets  are  monotonous  and  show  no  deviation  from  the  K=0  for  the  Ky,  implying  there  is  no  side  force.  This 
behavior  is  consistent  for  all  the  test  cases. 

On  the  other  hand,  weak  interaction  is  observed  when  the  jet  is  injected  into  the  leeward  side  (Bank=0).  At  high 
angles  of  attack,  the  Kz  curves  show  a  sudden  dip  in  the  range  of  14  to  18  degrees  of  AoA,  for  example,  at 
Mach=2.3  and  pressure  ratio  150.  A  notable  feature  is  that  the  Ky  curves  show  deviation  from  Ky=0  at  AoA  greater 
than  7  deg.  This  means  side  force  is  generated  even  though  the  flow  is  supposed  to  be  synnnetric.  This  prompted 
close  scrutiny  as  to  what  is  happening  for  the  leeward  jet  flow. 

In  Figs.  19  through  21,  jet  effectiveness  factor  for  normal  and  side  force  components,  Kz  and  Ky  respectively,  is 
displayed  for  Mach  2.3,  PR  150,  and  leeward  jet  which  is  the  same  case  shown  in  Fig.  6.  As  opposed  to  the 
windward  jet,  the  leeward  jet  interacts  with  the  body  vortex  when  the  angle  of  attack  reaches  a  certain  value.  It  is 
indeed  observed  from  the  schlieren  pictures  that  the  onset  of  body  vortex  travels  upstream  and  interacts  with  the 
bow  shock.  When  this  happens,  the  bow  and  barrel  shcok  suddenly  become  fuzzy  and  seem  to  vibrate,  causing  the 
flow  asymmetrical.  The  amount  of  side  force  is  not  negligible  and  is  as  much  as  20  percent  of  the  normal  force  at 
low  pressure  ratio. 


Figure  18.  Schlieren  graph  at  Mach  2.3  with 
leeward  jet. 


Figure  19.  Kz  as  a  function  of  AoA  at 
M=2.3  with  leeward  jet. 
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Figure  20.  Ky  as  a  function  of  AoA  at  21.  Ky  as  a  function  of  AoA  at 

M=2.3  with  leeward  jet.  M=2.3  with  leeward  jet. 


In  Figs  22  and  23,  side  force  factor,  Ky,  is  displayed  for  varied  Mach  numbers.  Fluctuations  are  measured  at  low 
Mach  number  with  leeward  jet  as  shown  Fig.  22.  This  side  jet  wind  tunnel  tests  show  the  complementary  aspect  of 
experiment.  While  CFD  displays  detailed  flow  behavior  of  transverse  jet  interactions,  the  seemingly  random  features 
of  extra  yaw  at  high  angles  of  attack  can  only  be  ascertained  by  the  testing.  This  shows  still  a  restricted  role  of  CFD 
in  the  course  of  missile  aerodynamic  design. 


Figure  22.  Ky  as  a  function  of  AoA  Figure  23.  Ky  as  a  function  of  AoA 

with  leeward  jet  for  varied  Mach  no.  with  windward  jet  for  varied  Mach  no. 


4.  Conclusions 

Four  flow  simulations  mentioned  here  exemplifies  the  usefulness  of  CFD  for  illustration  of  flow  behavior.  The 
complex  nature  of  shock/boundary  interaction,  transverse  jet  interaction  and  jet  impingement  are  well  captured 
numerically  and  their  accuracy  is  validated  through  wind  tunnel  testing  or  a  special  set  up.  As  was  shown,  the  high 
angle  of  attack  side  jet  flows  were  not  well  perceived  prior  to  the  wind  tunnel  test,  exemplifying  a  need  for 
computation  and  its  validation. 

Future  activities  at  the  ADD  aero  lab  include  multidisciplinary  optimization  in  an  effort  to  accelerating  the  design 
cycle. 
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Conclusion 


•  Aerodynamic  design  process  at  ADD  is  introduced. 

•  CFD  applications  for  cylinder  ogive,  boat  tail  and  spike, 
VLS  internal  flow,  side  jet  effect,  are  shown  along  with 
validation. 

•  Various  validation  techniques  are  devised  for  different 
flows. 

•  CFD  contributes  greatly  during  the  missile  development 
phase,  giving  insight  into  the  flow. 

•  Coupling  with  multidisciplinary  optimization  is  our  next 
direction. 
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